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Abstract 

The Large Area Telescope (LAT) on-board the Fermi Gamma-ray Space Telescope is a pair-conversion 
telescope designed to survey the gamma-ray sky from 20 McV to several hundreds of GeV. In this energy 
band there are no astronomical sources with sufficiently well known and sharp spectral features to allow an 
absolute calibration of the LAT energy scale. However, the geomagnetic cutoff in the cosmic ray electron- 
plus-positron (CRE) spectrum in low earth orbit does provide such a spectral feature. The energy and 
spectral shape of this cutoff can be calculated with the aid of a numerical code tracing charged particles 
in the Earth's magnetic field. By comparing the cutoff value with that measured by the LAT in different 
geomagnetic positions, we have obtained several calibration points between ^6 and <~13 GeV with an 
estimated uncertainty of ~2%. An energy calibration with such high accuracy reduces the systematic 
uncertainty in LAT measurements of, for example, the spectral cutoff in the emission from gamma ray 
pulsars. 

Keywords: Cosmic-Rays, geomagnetic cutoff, absolute energy scale, Fermi Large Area Telescope 



1. Introduction 

The Large Area Telescope (LAT) is the primary 
instrument on the Fermi Gamma-ray Space Tele- 

1 Corresponding author 



scope, which was launched on June 11, 2008. The 
LAT is a 4x4 modular array of identical tow- 
ers, each comprised of a tracker/converter (TKR) 
and a calorimeter (CAL), surrounded by an anti- 
coincidence detector (ACD). Each TKR module 
contains 18 layers of tungsten foil and paired silicon 
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strip detector x-y tracking planes. Each CAL mod- 
ule contains a hodoscopic array of 96 CsI(Tl) scin- 
tillating crystals 8.6 radiation lengths deep. High 
energy gamma rays predominantly convert in the 
TKR into electron-positron pairs, and the energy 
of the developing electromagnetic shower is mea- 
sured in the CAL. To cover the energy range of 
the instrument the individual CAL crystals must 
be able to accurately measure energy deposits be- 
tween 2 MeV and 60 GeV. To accomplish this, each 
CAL crystal can measure the deposited energy in 
four different overlapping energy ranges, namely 
2-100 MeV, 2 MeV-1 GeV, 30 MeV-7 GeV, and 
30 MeV-70 GeV. 

The process of converting the scintillation light 
measured in each individual CAL crystal to an en- 
ergy of the incident photon can be roughly divided 
into two phases. The first phase is the calculation 
of the energy deposition in each crystal (hereafter, 
the "crystal energy"), which relies on calibration 
constants derived from the measured signals from 
sea-level muons and on-orbit cosmic ray protons to- 
gether with constants derived from an on-orbit elec- 
tronic charge-injection system. The second phase is 
event reconstruction, which is comprised of a set of 
algorithms that account for energy deposited in the 
TKR, leakage of the shower out the back and sides 
of the CAL, and energy lost in gaps between the 
towers. We developed the energy reconstruction 
algorithms from a detailed instrument model and 
Monte Carlo simulation based on GEANT4 [l[ . 

To verify the fidelity of the instrument model and 
Monte Carlo simulation, we performed an extensive 
beam test campaign in 2006 at CERN and the GSI 
heavy ion accelerator laboratory. To eliminate the 
handling and shipping risks to the flight detector 
subsystems, and to reduce cost and schedule im- 
pacts, these tests were not done on the full LAT 
but on a Calibration Unit (CU), which consisted of 
two fully-populated flight-spare TKR modules and 
three fully-populated flight-spare CAL modules. A 
more detailed description of the beam test as well 
as the energy reconstruction algorithms will be pre- 
sented in section [21 

To perform an in-flight verification of the ab- 
solute energy scale of the LAT it is necessary to 
find an astrophysical source with a spectral fea- 
ture whose absolute energy peak and shape are well 
known. Because sources in the LAT energy band 
exhibit smooth spectral forms, it is not trivial to 
satisfy this requirement. Furthermore, the LAT 
is far more sensitive than any previous gamma- 



ray telescope, so it is not possible to use existing 
measurements of celestial sources to perform this 
calibration. The LAT has shown that gamma-ray 
pulsars typically exhibit spectral cutoffs in the 1- 
10 GeV range 2] and, for bright pulsars, the spec- 
tral cutoff energy can be measured with statisti- 
cal uncertainty approaching 1% [jj. If the system- 
atic uncertainty on the energy scale in this energy 
range can be reduced to the few-percent level, pre- 
cise measurements of spectral cutoff can address de- 
tailed questions of the emission physics and geome- 
try in pulsars [3j . In the orbit of the Fermi observa- 
tory (inclination of 25.6° and altitude of 565 km), 
the geomagnetic cutoff in the cosmic-ray electron 
plus positron (CRE) spectrum can serve as such a 
source in the range between ~6 and ~13 GeV. Elec- 
trons and positrons are a good calibration source 
because they interact in the LAT CAL in the same 
way that photons do: they both produce electro- 
magnetic showers. The energy measurement bias 
coming from the slight differences between CRE 
and photons interactions in the LAT is much less 
than 1% 0]. The capabilities of the LAT to mea- 
sure the CRE spectrum have been well demon- 
strated 0, [B[ . The cutoff rigiditjQ can be predicted 
by numerically tracing particle trajectories in the 
Earth's magnetic field, and the comparison of the 
predicted and measured values provides the oppor- 
tunity to perform this validation. 

In order to measure the deviation between the 
reconstructed electron spectrum and the calculated 
value as a function of energy, this analysis is per- 
formed in several Mcllwain L intervals. The Mcll- 
wain L-parameter 0] is a parameter describing a 
set of the Earth's magnetic field lines, in particular 
those which cross the Earth's magnetic equator at 
a number of Earth-radii equal to the L- value. Mag- 
netically equivalent positions (from the standpoint 
of the incoming charged particle) around the world 
will by definition have the same Mcllwain L values, 
therefore making this parameter particularly con- 
venient for characterizing cutoff rigidities \7\. The 
orbital inclination of the Fermi orbit fixes the range 
of Mcllwain L values accessible for this analysis to 
1.00-1.72. 



1 Rigidity is defined as the particle momentum divided by 
its charge. The cutoff we are measuring in the CRE spectrum 
is not the vertical cutoff but rather averaged over all angles. 
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2. LAT energy calibration 

The absolute energy scale of the LAT is defined 
by comparing the signals in the CAL crystals with 
the amount of energy a Monte Carlo simulation in- 
dicates should be deposited by on-orbit relativis- 
tic protons. The details of the energy calibration, 
referred to as the "proton inter-range calibration", 
are given in Q . To perform this calibration, we first 
correct the observed CAL signals for electronic non- 
linearities (as measured by an electronic charge- 
injection process) and for position-dependent scin- 
tillation response (as measured by a direct calibra- 
tion with sea-level muons and on-orbit with pro- 
tons). We then compare the corrected CAL sig- 
nals with the distribution of deposited energies pre- 
dicted by a GEANT4 simulation of the on-orbit 
spectrum of primary cosmic-ray protons passing 
through the LAT. This calibrates the highest-gain, 
lowest-energy range of the CAL readout. The re- 
maining gain ranges [1] of the CAL readout are cal- 
ibrated by enforcing that adjacent gain ranges give 
the same measured energy in the regions of energy 
space in which they overlap [8j . 

The Galactic cosmic-ray (GCR) element abun- 
dance peaks could in principle be used to calibrate 
the higher energy ranges of the LAT; however, the 
scintillation efficiency of heavy ions in CsI(Tl) is 
not identical to that of electromagnetic showers by 
an element-dependent factor that is not well known 
or measured. We measured the relative scintillation 
efficiencies of sub-relativistic ions in CsI(Tl) crys- 
tals in a series of beam tests [9| , but we do not have 
sufficient confidence that these efficiencies apply at 
relativistic energies. Thus, we are unable to re- 
late the observed signals from GCR element peaks 
to electromagnetic shower energy depositions accu- 
rately enough to calibrate the LAT energy scale. 

The geomagnetic cutoff in the CRE spectrum 
(measured in the energy range ~6 to ~13 GeV) 
is a good alternative source to calibrate the higher 
energy ranges of the CAL readout. At these en- 
ergies, the maximum energy per crystal is of the 
order of ~1 GeV, and will be read out in the higher 
energy ranges of the CAL. As a reference, 1 GeV of 
energy per crystal is equivalent to 100 times what a 
minimum ionizioning particle releases per crystal. 

The LAT CU was also calibrated using sea-level 
cosmic-ray muons, but the absolute energy scale 
was then cross-checked against calibrated beam 
lines at CERN, showing an energy deposit system- 
atically larger than expected. Such a direct calibra- 



tion relies on the knowledge of the incoming beam 
energy (determined with a 1% accuracy), the ge- 
ometry of the beam line and the CU detector, the 
Monte Carlo of the electromagnetic shower devel- 
opment within the detector, as well as on a good 
control of the environmental effects that might in- 
fluence the CU response (temperature, humidity, 
exceedingly high particle rates). 

The energy reconstruction of the LAT (as well as 
of the CU) is based on three different algorithms: 
a parametric correction based on the barycenter of 
the shower, a fit to the shower profile taking into ac- 
count the longitudinal and transverse development 
of the shower, and a maximum likelihood fit based 
on the correlation between the total deposited en- 
ergy, the energy deposited in the last layer of the 
CAL and the number of tracker hits. For each 
event, the best energy reconstruction method is se- 
lected by means of a classification tree analysis de- 
scribed in We would like to emphasize that 
at the energies considered in this analysis (i.e. be- 
tween ^6 and ~ 13 GeV) the correction factors for 
losses due to leakage out of the detector are reason- 
ably small (on average of the order of 30%) when 
compared to higher energies (of the order of 50% for 
100's of GeV) and therefore the method described 
here is testing a mixture of the crystal energy cali- 
bration and the leakage corrections. 

Following the beam test campaign, the GEANT4 
simulation of the CU was updated to reflect the best 
available description of the CU detector, the beam 
line, and the particle interactions. Most notably, 
the routines describ ing: the Landau-Pomeranchuk- 
Migdal effect (LPM 10]) in electromagnetic shower 
developments were updated after finding that ini- 
tial GEANT4 implementations were not providing 
a satisfactory description of the shower longitudinal 
development. The energy resolution for all three 
available energy reconstruction methods was mea- 
sured to be consistent with expectations. However 
despite these improvements, the reconstructed en- 
ergy in the CU was consistently higher than the 
beam energy by ~9%, on average, with further fluc- 
tuations of ~5%, depending on incoming beam en- 
ergy (5-282 GeV), angle (0-60 degrees) and posi- 
tion of the crystal within the shower 0, HH . 

The energy scale factor derived from the direct 
CU calibration was not applied to the LAT because 
of the differences in the beam test and on-orbit en- 
vironment, particularly the distribution and rate of 
the particles incident on the calorimeter, and the 
temperature variations for the CU during the data 
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taking, which were considered responsible for the 
different LAT and CU energy scales. 

3. Particle Tracing 

Geomagnetic cutoff rigidities can be obtained by 
tracing cosmic-ray trajectories in a model of the 
Earth's magnetic field. The standard mathemati- 
cal description of the Earth's magnetic field is given 
by the International Geomagnetic Reference Field 
models (IGRF) [l2j]. These models consist of the 
Gauss coefficients defining the spherical harmonic 
expansion of the geomagnetic potential up to a 
given order. In the IGRF-11 model, used for this 
analysis, the maximum multipole moment of the 
expansion is the 13 th . The IGRF model is com- 
puted by the participating members of the Interna- 
tional Association of Geomagnetism and Aeronomy 
(IAGA) Working Group V-MOD [lj] which is sup- 
ported by the organizations involved in operating 
magnetic survey satellites, observatories, magnetic 
survey programs and World Data Centers. IGRF- 
11 provides a definitive main field model for epoch 
2005, a main field model for 2010, and a linear pre- 
dictive secular variation model for 2010-2015. We 
use the particle trajectory tracing code (hereafter 
tracer) developed by Smart and Shea 0] and the 
IGRF-11 in this analysis to measure the geomag- 
netic cutoff. 

Since it is difficult to trace the trajectory of an 
incoming particle through the magnetic field and 
expect it to intersect the precise location desired, it 
is more efficient to calculate the trajectory in the re- 
verse direction. Therefore the starting point of the 
trajectory is given by the geographic coordinates, 
altitude and orientation of the spacecraft and the 
trajectories are propagated using IGRF-11 [12 1. 

At each of a grid of locations covered by the Fermi 
orbit, evaluated every 15 seconds from Aug 2008 
through July 2009, we simulated an ensemble of 
test particles. We selected energies of the test par- 
ticles according to the power-law spectrum mea- 
sured by the LAT [3] , and generated both electrons 
and positrons in the abundance ratio measured by 
PAMELA 13]. We use the tracer code to determine 
for each test particle whether it could have orig- 
inated from outside the geomagnctosphcrc. Tra- 
jectories that eventually intersect the Earth's at- 
mospher^l are rejected and labeled as forbidden 



Top view 




2 Distance taken to be 20 km from the Earth's surface. 



10.99 GV/ce- 
(vertical cutoff) 



Figure 1: Illustration of e — and e + of varying energies traced 
in the Earth's magnetic field from a point at longitude and 
latitude (0,0) and 565 km above the Earth's surface and 
viewed from the North Pole. The vertical cutoff is taken 
to be the lowest rigidity which the particle needs to escape 
when traced vertically through the Earth's magnetic field. 
As the rigidity of the particle decreases the amount of geo- 
magnetic bending increases and when the trajectory inter- 
sects the Earth it is taken to be forbidden (i.e of secondary 
origin). The cutoff values are different for positively and 
negatively charged particles because the Earth's magnetic 
field is not a perfect dipole. Particles labeled as escaped are 
taken to be of Galactic origin. Earth and trajectories are to 
scale. 



(secondaries hereafter). Trajectories that reach 20 
Earth radii are accepted as escaped (or of Galac- 
tic origin, primaries hereafter) . Figure [T] illustrates 
a few sample trajectories, both forbidden and al- 
lowed. 

We selected electron and positron data from the 
first year of Fermi LAT data taking (Aug 2008 
through July 2009) using the methods described in 
0] for the 100 MeV to 100 GeV energy range. In 
figure [H are the distributions of the reconstructed 
angle with respect to local zenith (upper panel) and 
azimuth (lower panel) for the flight data (black) and 
for the tracer output (red). The overall agreement 
between flight data and the tracer output is very 
good, demonstrating that the angular distributions 
have been well described. 



4. Analysis 

4-1- Estimating the fraction of primary Cosmic 
Ray electrons 

The population of CREs in low Earth orbit is 
a mixture of primary and secondary cosmic rays, 
where primary implies of extraterrestrial origin 
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Figure 3: Azimuthal distribution for flight data, tracer and secondaries in the Mcllwain L interval 1.0 < L < 1.14. The black 
distribution is the linear combination of tracer and secondaries which best fits the data. Each panel depicts the distributions in a 
given energy interval (labeled in the panel). In the upper left panel, where the secondary population is dominant, a maximum is 
evident in the westward direction where the secondary positrons are dominant. With increasing energy, the primary population 
gradually increases in the eastward direction. The LAT acceptance and rocking profile have been convolved, which in turn 
effects the overall shape of the azimuth distribution. The cutoff rigidity in this Mcllwain L range is ~13 GeV. The fraction of 
primaries evaluated from this template fitting is shown in figure [4] 



and secondary implies a mixture of splash and re- 
entrant particles produced in the interaction of pri- 
maries in the Earth's atmosphere. The trajectories 
of the secondary particles are extremely difficult to 
simulate reliably. By definition our tracing code 
only provides information on the escaped (i.e. pri- 
mary) particles. As a consequence it is necessary 
to estimate the fraction of secondary particles, as 
a function of energy, from the flight data. This 
constitutes one of the most delicate aspects of the 
analysis because the cutoff rigidity is found by fit- 
ting the spectrum (as will be described in detail 
in section 14.21) and its shape is influenced by the 
fraction of secondaries. 

The azimuthal distribution in Earth centered co- 
ordinates of the secondary population is different 
from that of the primary particles. We can exploit 
this fact and perform a template fit to identify the 



fraction of each population. To perform this task 
it is first necessary to choose the appropriate tem- 
plates to describe the populations. It is safe to as- 
sume that the population at low energy in flight 
data (E<<E C , where E c is the energy corresponding 
to the cutoff rigidity) is predominately composed of 
secondary particles, and can therefore be used as a 
template. The output from the tracer code provides 
the template for the primaries. Figure [3] shows an 
example of the template fitting performed to esti- 
mate the fraction of primaries in the interval 1.0 < 
Mcllwain L < 1.14. The gradual transition from a 
pure sample of secondary particles (top left panel, 
for energies between 2 and 4 GeV) to that of only 
primaries (bottom right panel, for energies between 
16 and 18 GeV) is evident. For reference, the cut- 
off energy in this Mcllwain L interval is ~13 GeV. 
For simplicity, only four energy intervals are shown. 
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Figure 2: The reconstructed angle with respect to local 
zenith (upper panel) and azimuth (lower panel) for data and 
tracer. In the lower panel, north is at 0° , east at 90° , south 
at 180° and west at 270°. The LAT acceptance for elec- 
trons and positrons as well as the rocking profile have been 
convolved in these distributions and influence their shape. 
Both are averaged over the orbit and for energies greater 
than 30 GeV in order to compare the distribution of pri- 
mary CREs. There is an overall good agreement between 
data and tracer. 

However the analysis has been performed in 16 over- 
lapping energy intervals. The resulting value for the 
fraction of primaries as a function of energy for this 
same Mcllwain L interval is shown in figure 21 

4-2. Measuring the cutoff rigidity 

The primary eT + e + spectrum can be parame- 
terized by 

dN(E) cE~ r 
dE (1 + (E/E c )~ 6 ) [ ' 

Where T is the spectral index and E c is the energy 
corresponding to the cutoff rigidity. We fit the pri- 
mary CRE spectrum (both data and tracer) with 
equation Q] to get the value of the cutoff rigidity. 



Figure 4: Fraction of primary electrons and positrons 
as a function of energy for the Mcllwain L interval 
1.00 < L < 1.14 

To reconstruct the CRE primary spectrum it is 
first necessary to remove the secondary population 
from the count spectrum and then correct for the 
hadronic contamination. The residual contamina- 
tion from hadrons, h, is estimated by applying the 
CRE event selection to the on-orbit simulation used 
by the Fermi LAT collaboration. We found that 
this residual hadron contamination ranges between 
^15% below the cutoff and ~5% above with an es- 
timated uncertainty on the absolute hadron flux of 
20%. The on-orbit simulation is based on data from 
CR experiments and includes all the components of 
the Galactic cosmic rays as well as the re-entrant 
and splash Earth albedo particles. A detailed de- 
scription of this simulation can be found in Q • The 
residual contamination due to secondaries, s, is es- 
timated as described in section BTTl The count spec- 
trum is multiplied by the purity factor defined as 
P.F = (1 — h) ■ (1 — s) in each energy bin. This 
background-subtracted count spectrum is divided 
by the effective geometry factor (EGF) [4| and the 
width of the energy interval, to obtain the final 
spectrum. 

The tracer count spectra are reconstructed by re- 
quiring the escaped condition and dividing by the 
width of the energy interval. To properly compare 
the cutoff rigidities, the tracer spectra have been 
binned in measured energy (i.e., they have been 
convolved using the energy resolution of the Fermi 
LAT for electrons and positrons |4j). No efficiency 
correction is required for the tracer spectra. To 
measure the spectra in Mcllwain L intervals an ex- 
tra cut specifying the interval is applied to both 
data and tracer. 
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Figure 5: Energy spectrum of data and tracer for the Mcll- 
wain L interval 1.0<L<1.14. The black line is a fit to the 
data using equation [l] and the resulting cutoff energy has a 
value of 13.27±0.10 GeV, and is indicated by the dot dashed 
line in the figure. The lower panel depicts the ratio of the 
two spectra. 

A comparison between the counts spectra mea- 
sured from flight data and from the tracer output 
is shown in figure [5j The shapes of the spectra are 
in good agreement, indicated by the ratio of the two 
depicted in the lower panel of the figure. Both the 
reconstruction and fitting procedures are applied to 
each Mcllwain L interval considered in this analysis. 

4-.S. Assessment of the systematic uncertainties 

The main sources of systematic uncertainty in 
this analysis are: the choice of the energy interval in 
which to fit the spectrum, the choice of the energy 
interval for the secondary template, the accuracy 
of the geomagnetic field model, and the accuracy of 
the Monte Carlo simulations used for the estimation 
of the residual hadron contamination as well as the 
evaluation of the EGF. 

To estimate the systematic uncertainty from the 
choice of energy interval used for the spectral fits, 
we performed a series of fits over a range of energy 
intervals narrower and broader than the chosen in- 
terval and calculated the root-mean square (rms) 
of the distribution of resulting values of the cutoff 
energy. The rms was found to be no greater than 
0.8% for all the Mcllwain L intervals. 



Since the key ingredients in the calculation of the 
fraction of CRE primaries are the secondary and 
primary templates, it is crucial to validate them. 
The tracer output (used as the template for the 
primaries) can be validated directly with the Fermi 
LAT flight data, as was described in section [3J 
The secondary template is more difficult to validate 
given that the output from the tracer code does not 
provide any information regarding the secondary 
populations. As stated in section I4.ll it is safe to 
assume that the population at low energy in flight 
data (E<<E C ) is predominately composed of sec- 
ondary particles. However the choice of the energy 
interval is not well determined, and can therefore be 
a source of uncertainty. This choice can in fact in- 
fluence the shape of the spectrum and thus the final 
value for the cutoff rigidity. To investigate its effect, 
we chose several different energies ranges in which 
to define our secondary templates (for each Mcll- 
wain L interval) by varying both the width of the 
energy interval as well as the distance from the cut- 
off energy. With the resulting fraction of primaries 
obtained from each template we reconstructed the 
spectrum and obtained the cutoff rigidity following 
the procedure described in section l4~2l The spread 
of the ratio of the cutoffs (found to be no greater 
than ~ 2% for all Mcllwain L intervals) gives an 
estimate for the uncertainty due to the secondary 
template. 

The accuracy of the cutoffs derived from the par- 
ticle tracing code is limited by the uncertainty of 
the geomagnetic field model. Several cross checks 
on the accuracy of the predicted rigidity cutoffs 
have been performed using satellite experiments, 
in particular the cosmic-ray isotope experiment 
HEAO-3 C2 14] (C2) and the Mass Spectrom- 
eter Telescope (MAST) on the Solar Anomalous 
and Magnetic Particle Explorer (SAMPEX) space- 
craft [ljf. The HEAO-3 satellite flew between 1979 
and 1981 at an altitude of 496 km and inclination 
of 43.6°. From the comparison between the exper- 
imental and computed cutoffs of oxygen nuclei in 
the 5 GeV/n range it was found that the calcu- 
lated cutoffs were systematically ^3-5 ± 2% higher 
than the measured ones [l4{. The SAMPEX mis- 
sion was launched into an orbit with 82° inclination 
with an average altitude of ^600 km and operated 
from 1992 to mid 2004. They performed a com- 
parison between measured cutoff rigidities of 0.3 
to 1.7 GeV protons during geomagnetically quiet 
times with those predicted by the Smart and Shea 
tracer code (based on IGRF-10) and found that the 
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latter systematically overestimate the data by 



14% [16 1 



Directly comparing these findings with the mea- 
surements presented here is difficult because in both 
the SAMPEX and HEAO-3 cases the measured cut- 
off rigidities are much lower than those appropriate 
for Fermi. Smart and Shea in 17[ assert that the 
accuracy of the model improves in regions closer 
to the geomagnetic equator (i.e. for larger cutoff 
rigidities). According to their analysis, it is reason- 
able to assume that the bias in the tracer output 
used in this analysis is no greater than ~3-5% for 
the Mcllwain L regions spanned by the Fermi ob- 
servatory. However, this is just an extrapolation 
based on the measurements performed by HEAO-3 
and SAMPEX at lower rigidities and a definitive 
answer on the value of this bias at the Fermi rigidi- 
ties remains unknown. 

The uncertainty of the hadronic component of 
the Monte Carlo simulations used to evaluate the 
residual contamination is estimated to be no greater 
than 20%. Its contribution to the overall systematic 
uncertainty in the final measurement of the cutoff 
is small (0.5%) because the contamination is always 
below 20% in every Mcllwain L interval. The as- 
sessment of the systematic uncertainty due to the 
EGF is described in detail in Q and for this energy 
range has a value of 10-15%. We studied its ef- 
fect on the final value of the cutoff via a bracketing 
method. To do this we generated a set of EGFs that 
have a maximal impact on the final shape of the 
spectrum, reconstructed the spectrum using these 
new EGFs and analyzed the final value of the cutoff 
energy obtained with these new EGFs. The result- 
ing values varied by ±1%. 

4- 4- Cross check of method 

To demonstrate that this analysis has the sensi- 
tivity to measure deviations of several percent in 
the energy scale of the LAT, we repeated the anal- 
ysis on a test sample of traced particles with an 
added +5% shift in energy. The resulting measure- 
ment of the cutoff energy was 5.4± 0.2% higher, 
which is consistent with the bias we introduced 
in the input. Thus we conclude that indeed this 
method is capable of measuring energy scale errors 
of the magnitude we expect. 

Pre-launch tests of crystal boule samples 
showed [HI that the light yield of the CsI(Tl) crys- 
tals in the LAT calorimeter could be expected to 
decrease by ~1% per year from radiation damage 



Drift (%/yr) 


B 


1.49±0.05 


C 


1.60±0.04 


N 


1.50±0.06 


O 


1.46±0.03 


CRE 


1.90±0.90 



Table 1: Drift per year for four abundant GRC element 
(B,C,N, and O) peaks and the CRE cutoff energy measured 
in Fermi LAT flight data over the first year of operations. 



in the charged-particle environment of low Earth or- 
bit, primarily from trapped particles in the South 
Atlantic Anomaly. 

By analyzing the path length-corrected^] crystal 
energies of four abundant GCR elements, namely 
boron, carbon, nitrogen and oxygen, over time it is 
possible to verify this prediction with flight data. In 
fact, all four of these peaks show a similar linear de- 
crease per year. These results are consistent with 
the pre-launch predictions. An example of these 
peaks is shown in figure [5] As a further cross check 
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Figure 6: The path length-corrected crystal energies of the 
B,C,N, and O peaks corresponding to the first four months 
of Fermi LAT operation. The large apparent abundance of 
B (compared to secondary and primary Galactic cosmic-ray 
abundance) is due to secondary B produced by primary C 
interacting in the ACD and TKR. 

of the sensitivity of the method used in this anal- 
ysis, we have compared the cutoff rigidity values 



3 The energy deposit is normalized by cos 8, where is the 
incidence angle of the particle. 
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measured in the first 60 days to the last 60 days 
of the first year of operations. We found that the 
measured cutoff energy decreased by 1.9±0.9% over 
this time interval, which is consistent with the de- 
crease from radiation damage measured with the 
GCR element peaks, as can be seen in table [T] We 
note also that this provides further confirmation of 
the sensitivity of our method. 

5. Results 

We applied this analysis to a data set spanning 
the first year of LAT operations. We found that 
the measured cutoff energy exceeded the predicted 
cutoff energy by 1.026±0.005 (stat) ±0.025 (sys) 
in the 6 GeV to 13 GeV range. The systematic 
error is the sum in quadrature of our estimates in 
Section 231 excluding those from the uncertainty in 
the magnetic field model and particle tracing code. 
As shown in Figure [71 the ratio between measured 
cutoff and predicted cutoff is constant within the 
errors over this energy range. 
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Figure 7: Distribution of cutoff energy ratios as a function 
of Mcllwain L. The ratios are constant within the errors 
over the entire range. The quoted value for the ratio is the 
weighted mean. The statistical errors are represented by the 
black error bars, and the systematic errors are in gray. As a 
reference, the value for the cutoff energy in each Mcllwain L 
interval is also shown. 



6. Conclusions 

We provide a method for calibrating the LAT en- 
ergy scale based on the measurement of the geomag- 
netic cutoff in the primary CRE spectrum. The re- 
sults shown here confirm that the LAT calorimeter 
calibration procedure, based on measured electronic 



gains and non-linearities, observed proton signals, 
and enforced consistency between readout energy 
ranges, is accurate to the ~2% level. It is worth 
noting that the method developed here is the only 
such calibration which relies uniquely on flight data 
and uses electromagnetic showers in the LAT. Ac- 
cordingly, this is a very important cross check on 
other methods such as beam tests and calibrations 
with cosmic-ray heavy ions. 

Our conclusion relics, however, on the fidelity of 
the current IGRF and particle trajectory tracing 
code. If we take the measurements of HEAO-3 C2 
and SAMPEX to indicate a bias of the order of +3- 
5%, the absolute energy scale of the LAT may have 
to be decreased by 5-7%. This conclusion, if con- 
firmed, would make the energy scale more consis- 
tent with the beam test results and the Monte Carlo 
simulations of the CU. This would suggest that the 
on-orbit energy calibration technique for the LAT 
should be revised. However, the cross checks re- 
ported in 14 , [l6[ on the predictions of the tracer 
code were both performed where the cutoff rigidi- 
ties are much lower than those measured by the 
LAT and in geomagnetic regions where the trajec- 
tory tracing calculations are less accurate 0- It is 
therefore unclear whether this bias in the theoret- 
ical cutoffs also applies to the results presented in 
this work. We therefore prefer to determine the en- 
ergy scale of the LAT using the results of the proton 
inter-range calibration, which are confirmed by the 
independent measurements discussed in this paper. 
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